Further  Evaluation  of  Methods  for 
Producing  Desired  Colors  on 
CRT  Monitors 


ASC  99-0581 


David  L.  Post,1*  Christopher  S.  Calhoun2 

1  Air  Force  Research  Laboratory.  Wnght-PaUerson  Air  Force  Base.  Ohio  45433  7022 
3  Logicon  Technical  Servicer.  Inc..  Dayton.  Ohio  45437-7258 


Received  13  March  IW;  accepted  16  July  1WV 


Abstract:  Three  experiments  were  performed  to  evaluate 
methods  for  predicting  the  luminances  and  chromaticity 
coordinates  produced  by  color  CRT  monitors,  given  known 
inputs.  Linear  and  logarithmic  versions  of  PLCC  and 
PLVC.  plus  Be  ms.  Motto,  and  Gorzvnski's  power  Junction, 
were  tested.  Estimates  are  provided  for  the  number  of  CRT 
measurement  points  needed  to  maximize  each  method's 
predictive  accuracy.  Correcting  for  unintended  light  from 
the  monitor  is  shown  to  improve  accuracy  substantially  for 
a  case  involving  a  seemingly  small  amount  of  light.  Rems  cl 
al.'s  characterization  technique,  which  involves  measuring 
the  monitor’s  neutral  [toini.  is  shown  to  yield  the  same 
accuracy  as  conventional  characterization  while  reducing 
the  number  of  measurements  required,  and  to  yield  im¬ 
proved  accuracy  when  correction  for  uninteiuled  light  is 
needed  but  impractical.  The  accuracies  of  the  predictive 
melliotls  are  conqtared  and  recommendations  for  their  use 
are  provided.  O  201)0  John  Wiley  &  Sons.  Inc.*  Col  Res  Appl.  25. 
00-104.  2000 

Key  words:  characterization:  color;  colorimetry:  CRT:  dis¬ 
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INTRODUCTION 

In  an  earlier  article.'  we  compared  several  methods  for 
predicting  the  colorimetric  behavior  of  digital-graphics  sys¬ 
tems  that  display  output  on  color  CRT  monitors  (DG/CRT 
systems).  We  found  that  two  methods  we  called  piecewise 
linear  interpolation  assuming  constant  chromaticity  coor¬ 
dinates  (PLCC)  and  piecewise  linear  interpolation  assum¬ 
ing  variable  chromaticity  coordinates  (PI .VC)  gave  the 
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most  accurate  predictions.  PLCC  uses  the  assumption  that 
the  chromaticity  coordinates  of  the  DG/CRT  system's  red. 
green,  and  blue  color  channels  (R.  G.  and  R)  arc  constant — a 
simplification  that  makes  PLCC  easy  to  implement  in  soft¬ 
ware.  PLVC  omits  the  constant  channcl-chromaticity  as¬ 
sumption  and  yields  more  accurate  predictions  of  chroma¬ 
ticity  at  the  expense  of  programming  ease.  (PLVC  is 
described  in  more  detail  in  Appendix  A.)  We  also  found  that 
16  measurements  of  each  color  channel,  spaced  equally 
between  the  smallest  digital-to-analog  converter  (DAC) 
value  that  can  be  measured  reliably  and  the  largest  possible 
DAC  value,  were  sufficient  to  maximize  PLCC  and  PL  VC’s 
predictive  accuracies. 

More  recently,  other  researchers  have  reported  new  ideas 
that  deserve  investigation.  Por  example.  Lucusscn  and  Wal- 
raven-  have  demonstrated  a  variation  on  PLCC  in  which  log 
luminance  is  interpolated  as  a  function  of  log  DAC-value — 
a  method  we  label  PLGCC.  This  refinement  might  be  ex¬ 
pected  to  yield  better  accuracy  than  PLCC.  because  a  log- 
log  transform  tends  to  linearize  DAC-value  — *  luminance 
transfer  functions  (the  degree  of  success  depends  on  the 
transfer  functions').  Applying  the  same  modification  to 
PLVC.  yielding  a  method  we  label  PI.GVC,  might  improve 
its  accuracy,  too. 

Another  important  development  has  been  Berns.  Motta, 
and  Gor/ynskiV  introduction  of  a  predictive  method  that 
uses  a  power  function  and  has  become  a  de  facto  CIE 
recommendation.4  Berns'  later  named  this  the  gain-offset- 
gamma  (GOG)  model.  Berns  el  al.’s  analysis  indicated  that 
GOG  is  as  accurate  as  PLCC,  but  requires  a  total  of  only 
eight  measurements  to  achieve  this  accuracy,  when  it  is  used 
in  conjunction  with  a  new  characterization  technique.  This 
technique  involves  measuring  the  tristimulus  values  of  the 
DG/CRT  system's  neutral  point  (N)  at  five  DAC  triplets  for 
which  l)s  -  l)(;  -  D„  T  0  (where  PK  represents  an  R  DAC 
value,  etc.)  and  performing  one  measurement  each  of  R,  G, 
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and  B  alone  to  determine  the  channel  chromaticities.  Given 
the  assumption  of  constant  channel  chromaticity,  the  results 
from  the  individual  R.  G.  and  B  measurements  can  be  used 
to  decompose  the  five  measurements  of  N  into  their  as¬ 
sumed  constituent  RGB  luminances,  yielding  five  lumi¬ 
nance  values  for  all  three  channels  with  only  five  measure¬ 
ments.  This  procedure  requires  fewer  measurements  than 
the  more  common  characterization  technique,  which  in¬ 
volves  measuring  R,  G.  and  B  alone  at  multiple  DAC  values 
with  the  other  two  channels  set  to  zero,  and  is  easy  to 
automate  in  software. 

Still  another  important  innovation  has  resulted  from  the 
realization  that  measuring  a  given  screen  location  may  yield 
a  nonzero  reading,  even  though  the  DAC  values  associated 
with  that  location  are  zero,  the  measuring  environment  is 
dark,  and  the  measuring  instrument  is  zeroed.  This  situation 
can  produce  large  prediction  errors  and  has  several  possible 
caascs:  (I)  the  monitor’s  black  level,  or  DC  offset,  is  too 
high:  (2)  the  cutoff  voltages  produced  by  the  graphics  card’s 
DACs  are  too  high;  (3)  light  from  an  adjacent  screen  area 
backscattcrs  from  the  faceplate  onto  the  location  being 
measured  (a  phenomenon  known  as  internal  or  interreflec- 
lion  flare);  (4)  electrons  from  an  adjacent  screen  area  back- 
scatter  from  the  shadow  mask  or  aluminized  screen  and 
excite  the  measuring  location:  and  (5)  light  from  a  source 
outside  the  intended  measuring  location  is  collected  by  the 
measuring  instrument's  optics  (a  phenomenon  called  instru¬ 
ment  flare).  Howard'-7  was  apparently  the  first  author  to 
correct  for  such  problems,  which  we  refer  to  collectively  as 
black  light,  and  others  have  followed  suit.4-*-*-*  Del  Barco, 
Diaz,  Jimdncz,  and  Rubifio.1"  in  particular,  have  provided  a 
detailed  analysis  that  shows  the  improvement  in  prediction 
accuracy  that  black  correction  can  produce. 

The  present  experiments  were  designed  to  update  our 
earlier  evaluation'  by  trying  out  the  ideas  discussed  above. 
Specifically,  we  sought  to  quantify  the  effects  of  black 
correction  and  neutral-point  characterization  on  the  accura¬ 
cies  of  the  aforementioned  predictive  methods,  obtain  re¬ 
fined  estimates  of  the  optimal  number  of  characterization 
measurements,  and  compare  results  across  methods.  Our 
overall  purpose  is  to  provide  information  that  is  useful  in 
selecting  a  method,  obtaining  the  best  performance  from  it, 
and  developing  more  accurate  methods. 


APPARATUS 

Video  Generation 

Video  signals  were  produced  by  a  Univision  model 
UDC7032A-I0C  graphics  board  operating  under  the  control 
of  our  software.  Each  image  plane  on  the  graphics  board 
generates  1024  X  1024  8-bit  pixels  at  a  60-Hz  noninterlaced 
refresh  rate.  The  board  is  equipped  with  an  Analog  Devices 
model  ADV7I50  135  MHz  video  RAM-DAC,  which  con¬ 
tains  a  3  X  256  X  10-bit  lookup  table  that  addresses  three 
10-bit  DACs.  Thus,  the  largest  available  DAC  value  for 
each  color  channel  (/),ml>)  is  1023. 


Colorimetric  Instrument 

Colorimetric  measurements  were  performed  with  an 
LMT  model  1210  filter  colorimeter,  equipped  w  ith  an  RS- 
232  interface  to  permit  external  computer  control.  The 
LMT's  optical  head  contains  no  lenses  for  focusing  light 
and.  therefore,  docs  not  yield  luminance  measurements; 
instead,  it  gathers  light  using  a  diffusing  integrating  surface 
and  measures  illuminance.  Therefore,  when  measuring  with 
the  LMT.  we  placed  its  optical  head  flush  with  the  CRT 
monitor’s  faceplate  to  maximize  the  signal.  The  LMT  lu¬ 
minance  measurements  we  report  subsequently  were  ob¬ 
tained  using  an  illuminancc-lo-luminancc  conversion  con¬ 
stant,  computed  by  measuring  the  display  at  DH-DG  =  Dn 
-  1023  with  the  LMT  and  a  spcctroradiomctcr  and  ratioing 
the  measurements. 

After  loading  each  new  set  of  DAC  values  into  our 
graphics  board,  we  allowed  90  s  to  elapse  before  commenc¬ 
ing  measurement  to  reduce  the  effects  of  short-term  monitor 
instability.  Our  software  then  commanded  the  LMT  to  per¬ 
form  four  measurements  in  succession  and  return  the  aver¬ 
aged  results.  All  measurcmenis  were  performed  in  a  dark 
room. 

Display 

We  tested  a  19-inch  (48-cm)  Aydin  model  9000  color 
CRT  monitor.  It  was  adjusted  to  provide  a  1:1  aspect  ratio, 
in  keeping  with  the  graphics  board's  1024  X  1024  addres¬ 
sability.  and  to  yield  a  D(>*  peak  white  having  a 

luminance  of  approximately  72  cd/nr  at  f)H  —  Da  -  l)„  - 
l),^,  and  a  perceptual  black  when  the  entire  screen  was  set 
to  Dh  ~  Dc  =  Dn  =  0  and  observed  for  5  min  in  the  dark 
room.  The  monitor  was  warmed  up  24  h  prior  to  data 
collection  to  improve  its  temporal  stability. 

Measuring  Target 

The  measuring  target  was  a  solid  circle  drawn  at  the 
center  of  the  CRT  monitor's  screen  against  a  gray  back¬ 
ground  having  DAC  values  =  (590,  590,  590).  These  values 
produced  a  luminance  that  was  roughly  20'S  of  IV**,,.  Die 
circle  was  400  pixels  in  diameter  and.  therefore,  occupied 
approximately  12%  of  the  display’s  active  area.  This  mea¬ 
suring  target  and  background  arc  similar  with  those  used  by 
Bems  et  at.'  Before  each  measuring  sequence,  the  LMT’s 
optical  head  was  positioned  flush  with  the  faceplate  and 
centered  within  the  circle. 

PROCEDURE 

Our  data  collection  was  designed  to  produce  results  that  can 
be  compared  with  those  reported  by  Bems  et  a!.y  Our  basic 
procedure  was  to  perform  a  scries  of  characterization  mea¬ 
surements,  followed  immediately  by  measurement  of  a  set 
of  test  colors.  This  process  was  repeated  five  times.  The 
characterization  data  from  each  replication  were  used  to 
predict  the  luminances  and  chromaticity  eexirdinates  of  the 
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associated  test  colors,  using  various  methods,  and  two  de¬ 
pendent  measures  were  computed  for  each  test  color:  ( 1  > 
absolute  value  of  the  percent  luminance  error:  and  (2) 
chromaticity  error,  expressed  as  distance  on  the  CIF.  1976 
uniform  chroinaticity-scak  (UCS)  diagram. 

Characterization 

For  each  replication,  the  DG/CRT  system  was  character¬ 
ized  by  measuring  each  color  channel's  luminance  and  C1E 
1976  u’v’-chromaticity  coordinates  at  several  DAC  values 
with  the  other  two  channels  set  to  zero.  Specifically,  we 
measured  R.  (7,  and  R  at  49  DAC  values  each  that  were 
spaced  equally  between  the  smallest  DAC  value  that  could 
be  measured  reliably  (380,  specifically)  and  D max.  The  ad¬ 
vantage  of  measuring  49  points  in  this  way  is  that  one  can 
form  subsets  containing  25,  17.  13,  9,  7.  5.  4.  3.  or  2 
measurements,  all  of  which  arc  spaced  equally  between  and 
include  the  smallest  measurable  and  largest  possible  DAC 
values. 

We  also  measured  the  neutral  points  (/V)  associated  with 
the  aforementioned  DAC  values.  That  is.  we  measured  the 
49  DAC-value  triplets  for  which  DR,  =  DCI  =  DB  r  where 
i  denotes  one  of  the  49  DAC  values.  Each  N  measurement 
was  then  decomposed  into  its  constituent  RGB  luminances, 
yielding  a  second  characterization  dataset. 

Test  Colors 

Our  test  dataset  consisted  of  the  same  DAC  triplets  as 
Berns  el  ul .\  except  for  (0.  0.  0).*  Specifically,  we  mea¬ 
sured  the  DAC  values  0.  559.  755.  903.  and  1023  in  all 
possible  three-way  combinations,  except  for  (0. 0. 0).  yield¬ 
ing  124  measurements.  The  resulting  luminances  and  chro¬ 
maticity  coordinates  were  different  from  Berns  el  al.’s  pre¬ 
sumably.  because  our  monitor's  peak  luminance  (72  cd/nr’) 
was  somewhat  higher  than  theirs  (65  cd/m:)  and  because  of 
probable  differences  in  the  monitors'  phosphors  and  volt- 
age-to-luminance  transfer  functions. 

Predictive  Methods 

We  replicated  Berns  el  al.’s'  comparison  of  GOG  using 
the  conventional  approach  to  characterization  vs.  the  use  of 
neutral-point  measurements.  (We  made  a  correction  to  their 
neutral-point  characterization  procedure,  though:  it  is  dis¬ 
cussed  in  Appendix  B.)  We  also  tested  Pl.CC.  PLGCC, 
PLVC.  and  PLGVC  using  conventional  and  neutral-point 
characterizations.  The  chromaticity  coordinates  obtained  at 
D,ml,  for  each  channel  were  used  to  represent  the  coordi¬ 
nates  of  the  RGB  channels  for  GOG,  PLCC.  and  PLGCC. 


•  We  omilleU  this  triplet  because  the  complete  data  collection  effon 
included  full  replications  using  a  MAG  model  MX  1 71  monitor,  as  well  as 
a  Photo  Research  model  PR-703A/PC  spectruradiumeter.  which  could  not 
measure  the  luminances  at  (0.  0.  0).  Tire  results  from  these  replications 
matched  those  for  the  Aydin  I.MT  combination  in  all  important  respects, 
and  have,  therefore,  been  omitted  for  the  sake  of  brevity. 


which  assume  constant  channel  chromaticity.  PLVC  and 
PLGVC  used  the  coordinates  that  were  measured  at  each 
DAC  value. 

The  idea  of  using  neutral-point  characterization  with 
PLVC  (or  PLGVC)  might  seem  nonsensical  because  the 
decomposition  process,  as  we  have  described  it.  relies  on 
the  assumption  that  the  RGB  channel  chromaticities  are 
constant.  Therefore,  a  characterization  dataset  produced  by 
this  process  contains  R.  G.  and  B  measurements  that  have 
constant  chromaticity — at  least  implicitly — and  using  such 
a  dataset  with  PLVC  yields  the  same  results  as  PLCC. 

It  is  possible,  though,  to  combine  information  from  con¬ 
ventional  and  neutral-point  characterizations  to  obtain  a 
decomposed  neutral-point  characterization  dataset  that  con¬ 
tains  information  about  measured  changes  in  channel  chro¬ 
maticity.  The  procedure  is:  (1)  Use  PLVC  and  the  results 
from  a  conventional  characterization  to  compute  predicted 
RGB  DAC  values  for  each  A'  measurement.  (2)  Use  PLVC 
and  the  same  conventional  characterization  to  predict  the 
luminances  and  chromaticity  coordinates  of  the  individual 
R.  G.  and  B  DAC  values  computed  in  step  1;  and  (3) 
Substitute  the  DAC  values  at  which  each  A'  measurement 
was  performed  for  those  computed  in  step  I .  For  example, 
an  Af  measurement  at  DR  =  Da  —  D„  =  511  might  yield 
computed  DAC  values  of  (509.  508.  512);  the  predicted 
luminances  and  chromaticity  coordinates  for  (509. 0.  0).  (0, 
508.  0).  and  (0.  0.  512)  would  then  be  computed  and  the 
associated  DAC  values  changed  to  (51 1.  0,  0),  (0,  511,  0), 
and  (0.  0.  511). 

Step  3  is  essential  because,  if  it  is  omitted,  the  resulting 
characterization  dataset  is  purely  a  transformation  of  the 
original,  conventional  characterization  and  any  differences 
in  the  resulting  predictions  arc  due  only  to  changing  the 
interpolation  endpoints  (plus  possible  differences  in  round¬ 
ing  and  truncation).  The  correct  procedure,  however,  yields 
a  dataset  in  which  channel  chromaticity  varies  with  DAC 
value,  in  accordance  with  the  measurements  from  the  con¬ 
ventional  characterization,  but  which  allows  PLVC  to  pre¬ 
dict  the  measurements  from  the  neutral-point  characteriza¬ 
tion  perfectly.  In  our  implementation  of  this  decomposition 
method,  we  used  conventional  characterization  datasets  that 
contained  the  same  number  of  measurements  as  the  neutral- 
point  characterization  dataset  that  was  undergoing  decom¬ 
position. 

As  in  our  earlier  article,  our  implementations  of  PLCC 
and  PLVC  used  the  simplifying  assumption  that  each  chan¬ 
nel's  luminance  is  zero  at  zero  DAC  value.  Our  experience 
has  been  that,  if  the  monitor's  black  level  is  adjusted  to  yield 
a  perception  of  black  for  a  dark-adapted  viewer,  this  as¬ 
sumption  yields  reasonably  accurate  results  and  avoids  the 
problems  inherent  in  trying  to  measure  extremely  low  (and 
possibly  zero)  lighl  levels.  For  PLGCC  and  PLGVC.  how¬ 
ever.  the  assumption  cannot  be  used,  because  zero  DAC 
value  and  luminance  do  not  exist  in  the  log-log  domain;  that 
is.  log(0)  is  undefined.  Our  solution  to  this  problem  was  to 
use  the  slope  defined  by  the  two  smallest  measurements  for 
a  given  channel  to  extrapolate  for  values  below  the  mea¬ 
sured  range,  when  necessary. 
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TABLE  I.  Mean  errors  from  Experiment  1. 
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1.42 

0.61 

1.17 

48 

1.18 

0.00141 

1.15 

2.46 

0.92 

1.90 

0.63 

0.00084 

0.78 

1.41 

0.61 

1.18 

PLCC 

7 

1.84 

0.00179 

1.48 

2.85 

1.19 

2.19 

0.87 

0.00131 

1.11 

2.01 

0.87 

1.66 

9 

1.55 

0.00152 

1.27 

2.51 

1.03 

1.95 

0.67 

0.00102 

0.88 

1.65 

0.69 

1.39 

13 

1  26 

0.00140 

1.15 

2.40 

0.93 

1.87 

0.62 

0.00085 

0.77 

1.43 

0.61 

1.23 

17 

1.25 

0.00141 

1.16 

2.35 

0.93 

1.83 

0.63 

0.00088 

0.81 

1.43 

0.64 

1.22 

25 

1.17 

0.00135 

1.11 

2.34 

0.89 

1.80 

0.63 

0.00080 

0.75 

1.38 

0.59 

1.13 

49 

1.15 

0.00133 

1.10 

2.34 

0.88 

1.77 

0.65 

0.00076 

0.73 

1.34 

0.58 

1.09 

PLGCC 

2 

4.54 

0.00401 

3.46 

6.54 

2.72 

5.28 

6.70 

0.00368 

3.68 

6.34 

2.91 

4.90 

3 

1.01 

0.00159 

1.21 

2.67 

098 

2.25 

1.90 

0.00118 

1.01 

2.00 

0.83 

1.64 

4 

0.78 

0.00121 

1.00 

2.46 

0.79 

1.90 

1.05 

0.00052 

0.63 

1.14 

0.51 

0.92 

5 

1.01 

0.00125 

1.0*1 

2.26 

0.83 

1.77 

0.72 

0.00069 

0.70 

1.27 

0.55 

1.03 

7 

1.05 

0.00126 

1.06 

2.24 

0.84 

1.70 

0.70 

0.00065 

0.66 

1.23 

0.52 

0.99 

9 

1.10 

0  00130 

1.07 

2.28 

0.86 

1.81 

0.64 

000075 

0.72 

1.32 

0.56 

1.07 

13 

1.11 

0  00131 

1.09 

229 

087 

1.75 

0.64 

0.00073 

0.70 

1.32 

0.56 

1.05 

17 

1.15 

0.00133 

1.11 

226 

0.88 

1.77 

0.64 

0.00078 

0.74 

1.32 

0.59 

1.10 

25 

1.13 

0.00133 

1.10 

2.31 

0.88 

1.76 

0.64 

0.00077 

0.73 

1.34 

0.58 

1.10 

49 

1.13 

0.00132 

1.09 

2.32 

0.88 

1.74 

0.66 

000075 

0.72 

1.33 

0.57 

1.08 

PLVC 

7 

1  84 

000125 

1.06 

1.97 

0.85 

1.46 

0.80 

0.00104 

0.92 

1.76 

0.73 

1.43 

9 

1.55 

0.00097 

0.87 

1.66 

0.71 

1.19 

0.59 

0.00075 

0.71 

1.43 

0.56 

1.12 

13 

1.26 

0.00089 

0.80 

1.52 

0.64 

1.11 

0.56 

0.00059 

0.64 

1.21 

0.52 

1.00 

17 

1.25 

0.00087 

0.79 

1.50 

0.63 

1.10 

0.58 

0.00062 

0.67 

1.21 

0.54 

0.98 

25 

1.17 

000084 

0.77 

1  48 

0.62 

1.07 

0.58 

000056 

0.64 

1.20 

0.52 

0.99 

49 

1.15 

0.00082 

0.75 

1.47 

0.61 

1.10 

0.61 

0  00052 

0.63 

1.19 

0.51 

096 

PLGVC 

2 

4.54 

0.00380 

3.68 

6.30 

3.01 

5.19 

6.54 

0.00445 

4.31 

7.43 

3.56 

6.26 

3 

1.01 

0.00129 

1.15 

2.13 

0.95 

2.00 

1.85 

0.00143 

1.26 

2.80 

1.09 

2.63 

4 

0.78 

0.00065 

0.68 

1.41 

0.55 

1.14 

1.00 

0.00049 

0.66 

1.26 

0.57 

1.20 

5 

1.01 

0.00078 

0.74 

1.42 

0.59 

1.05 

0.68 

0.00048 

0.63 

1.14 

0.52 

1.02 

7 

1.05 

0.00075 

0.73 

1  42 

0.59 

1  08 

0.65 

0.00046 

060 

1.15 

0.49 

0.98 

9 

1.10 

0.00078 

0.73 

1.41 

0.59 

1.06 

0.59 

0.00052 

0.63 

1.21 

0.51 

1.00 

13 

1.11 

0.00080 

0.76 

1.44 

0.60 

1.07 

0.61 

0.00050 

0.61 

1.16 

0.50 

0.97 

17 

1.15 

0.00081 

075 

1  46 

0.61 

1.09 

0.60 

0.00054 

0.64 

1.18 

0.52 

099 

25 

1.13 

0.00082 

0.76 

1  47 

0.61 

1  08 

0.59 

0.00054 

0.64 

1.21 

0.52 

1.00 

49 

1.13 

0.00081 

0.75 

1.48 

0.60 

1.10 

0.61 

0.00054 

0.65 

1.18 

0.52 

0.99 

EXPERIMENT  I:  RESOEITION  AND 
CHARACTERIZATION  METHOD 

Our  first  experiment  was  a  partial  replication  of  Berns  et 
aU  extended  to  include  Pl.VC.  PI.GCC.  and  PLGVC.  Each 
method’s  predictions  of  the  124-color  test  dataset  were 
tested  using  two  repeated- measures  full-factorial  analyses 
of  variance  (ANOVAs)  in  which  the  error  tenn  for  each 
effect  was  the  interaction  of  that  effect  with  replications. 
One  ANOVA  used  absolute  percent  luminance  error  as  the 
dependent  measure,  and  the  oilier  used  u'v'  error.  In  both 
cases,  the  dependent  measures  were  averaged  across  the  124 
lest  colors.  The  independent  variables  in  each  ANOVA 
were  characterization  method  (conventional  vs.  neutral- 
point)  and  the  number  of  measurements  per  color  channel  in 
the  characterization  dataset  (which  we  tenn  resolution).  The 
significance  of  all  statistical  tests  was  assessed  using  an 
alpha  of  0.05. 

The  results  arc  summarized  in  Table  I  and  Fig.  I.  Mean 


values  of  and  AE*,,.  calculated  relative  to  W _ for 

each  replication,  have  been  included  in  Table  1  to  facilitate 
comparisons.  The  95th  percentiles  for  A E*uv  and  A E?lh  arc 
shown  also,  to  provide  information  about  the  error  ranges. 

gog 

For  GOG.  we  planned  originally  to  test  resolutions  2-49. 
For  resolution  2.  though,  we  had  problems  solving  for  the 
coefficients  used  in  the  GOG  power  functions.  Eventually, 
we  realized  that  this  w  as  because  of  an  inherent  problem — 
explained  in  Appendix  C — that  arises  wiien  the  character¬ 
ization  dataset  includes  a  measurement  at  D^.  Here,  we 
w-ill  say  simply  that  measurements  contribute  little  or 
nothing  to  solutions  for  the  GOG  coefficients.  We  circum¬ 
vented  this  problem  by  substituting  measurement  48  for 
measurement  49  in  our  GOG  tests.  As  a  result,  the  highest 
resolution  we  could  test  for  GOG  was  48. 
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FIG.  1.  Mean  |%Y  and  u'v'  error  from  Experiment  1  for 
each  predictive  method  as  a  function  of  resolution  and  char¬ 
acterization  method. 


Tlie  ANOVA  on  luminance  emir  showed  that  ihe  main 
effects  of  characterization  and  resolution  are  significant  and 
account  for  53.3  and  17.7%  of  the  variance,  respectively. 
Neutral-point  characterization  produced  less  error  than  con¬ 
ventional  characterization,  and  increasing  the  resolution 
from  2  to  3  produced  a  much  larger  error  reduction  than  any 
other  resolution  change. 

The  ANOVA  on  u'v'  emir  showed  that  the  main  effects 
of  characterization  and  resolution,  plus  their  interaction,  are 
significant  and  account  for  43.9.  23.1.  and  7.0%  of  the 
variance,  respectively.  Again,  neutral-point  characterization 
produced  less  emir  than  conventional  characterization,  and 
increasing  the  resolution  from  2  to  3  produced  the  largest 
error  reduction.  The  interaction  appears  to  reflect  the  fact 
that  the  initial  drop  is  much  larger  for  conventional  charac¬ 
terization. 


The  effects  of  resolution  are  fairly  consistent  but  quite 
modest  beyond  resolution  3.  Therefore,  we  conclude  that 
there  is  little  to  be  gained  by  using  resolutions  beyond  3 
with  GOG.  Also,  neutral-point  characterization  yields  less 
error  than  conventional  characterization. 


PI.CC 

For  Pl.Ce.  we  tested  resolutions  7-49.  because  our  pre¬ 
vious  work1  indicated  that  lower  resolutions  would  not  yield 
competitive  accuracy.  The  ANOVA  on  luminance  error 
showed  that  the  main  effects  of  characterization  and  reso¬ 
lution.  plus  their  interaction,  are  significant  and  account  for 
70.2.  16.0,  and  4.6%  of  the  variance,  respectively.  Neutral- 
point  characterization  produced  less  error  than  conventional 
characterization.  Increasing  resolution  caused  a  larger  and 
completely  monotonic  emir  reduction  for  conventional 
characterization,  although  the  improvements  beyond  reso¬ 
lution  13  arc  trivial.  For  neutral-point  characterization,  error 
asymptotes  somewhere  around  resolution  13. 

The  ANOVA  on  u'v'  error  showed  that  the  main  effects 
of  characterization  and  resolution  are  significant  and  ac¬ 
count  for  65.9  and  27.7%  of  the  vuriance,  respectively. 
Again,  neutral-point  characterization  produced  less  error 
than  conventional  characterization.  Increasing  resolution 
produced  nearly  monotonic  emir  reductions,  although  the 
changes  beyond  resolution  13  are  modest. 

The  effects  of  resolution  are  fairly  consistent  but  very 
small  beyond  resolution  13.  Therefore,  we  conclude  that 
there  is  little  to  be  gained  front  using  resolutions  beyond  1 3 
with  PI.CC.  Also,  neutral-point  characterization  is  clearly 
preferable. 


PLGCC 

For  PLGCC.  we  tested  resolutions  2-49.  Ihe  ANOVA 
on  luminance  error  showed  that  the  main  effect  of  resolution 
and  the  characterization  .t  resolution  interaction  are  signif¬ 
icant  and  account  for  91.3  and  7.9%  of  the  variance,  respec¬ 
tively.  Increasing  resolution  caused  a  larger  and  more  con¬ 
sistent  decrease  in  error  for  neutral-point  characterization. 
As  a  result,  conventional  characterization  produced  lower 
error  at  resolutions  below  5.  but  neutral-point  characteriza¬ 
tion  is  better  at  the  higher  resolutions.  Error  for  conven¬ 
tional  characterization  reaches  a  minimum  at  resolution  4 
and  then  increases,  stabilizing  somewhere  around  resolution 
9;  emir  for  neutral-point  characterization  also  stabilizes 
around  resolution  9. 

The  ANOVA  on  u’v’  error  showed  that  the  main  effects 
of  characterization  and  resolution,  plus  their  interaction,  are 
significant  and  account  for  9.0.  X9.9.  and  0.3%  of  the 
variance,  respectively.  Neutral-point  characterization  yielded 
less  error  than  conventional  characterization  in  all  cases. 
Error  reaches  a  minimum  at  resolution  4  and  then  increases 
slightly,  stabilizing  somewhere  around  resolution  9.  At  the 
lower  resolutions,  there  are  minor  differences  in  the  effect 
of  resolution  for  conventional  vs.  neutral-point  character- 
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ization;  these  differences  appear  to  be  the  source  of  the 
smull  interaction. 

Both  types  of  error  seem  to  stahili/e  around  resolution  9, 
so  we  conclude  that  there  is  no  assured  advantage  to  using 
higher  resolutions  with  PLGCC.  Examination  of  the 
CIELUV  and  CIELAB  values  and  their  95th  percentiles  in 
Table  I  shows  that  a  case  can  be  made  for  using  resolutions 
4  or  5.  This  case  is  not  as  uncertain  as  it  might  seem, 
because  we  have  obtained  comparable  results  from  tests 
involving  a  different  monitor  and  colorimetric  instrument. 
We  conclude,  therefore,  that  resolution  4  or  5  is  best,  but 
recommend  9  for  those  who  prefer  a  conservative  alterna¬ 
tive.  Based  on  the  CIELUV  and  CIELAB  values,  we  also 
conclude  that  neutral-point  characterization  is  superior  to 
conventional  characterization  at  resolutions  of  3  or  more. 

PLVC 

for  PLVC,  we  tested  resolutions  7-  49.  The  ANOVA  on 
luminance  error  showed  that  the  main  effects  of  character¬ 
ization  and  resolution,  plus  their  interaction,  are  significant 
and  account  for  73.3.  13.8.  and  4.6%  of  the  variance, 
respectively.  Examination  of  these  effects  leads  to  the  same 
conclusions  we  reported  above  lor  PLCC.  This  outcome  is 
not  surprising,  because  our  earlier  article1  showed  that 
PLCC  and  PLVC  produce  nearly  identical  luminance  errors. 

The  ANOVA  on  u  V  error  showed  that  the  main  effects 
of  characterization  and  resolution,  plus  their  interaction,  arc 
significant  and  account  for  35.6.  54.6.  and  0.7%  of  the 
variance,  respectively.  Neutral -point  characterization  pro¬ 
duced  less  error  than  conventional  characterization  in  all 
cases.  Error  decreases  monolonically  with  resolution  for 
conventional  characterization,  but  seems  to  stabilize  around 
resolution  13  for  neutral-point  characterization. 

The  effects  of  resolution  .ire  consistent  but  rather  small 
beyond  resolution  13  for  conventional  characterization;  they 
stabilize  around  resolution  13  for  neutral-point  character¬ 
ization.  Therefore,  we  conclude  that  there  is  little  to  be 
gained  from  using  resolutions  beyond  13  with  PLVC.  Also, 
neutral-point  characterization  is  preferable. 

PLGVC 

For  PLGVC.  we  tested  resolutions  2-49.  The  ANOVA 
on  luminance  error  showed  that  the  main  effects  of  resolu¬ 
tion  and  the  characterization  x  resolution  interaction  are 
significant  and  account  for  91.6  and  7.5%  of  the  variance, 
respectively.  Examination  of  the  interaction  leads  to  tire 
same  conclusions  we  reported  above  for  PLGCC;  thus,  it 
seems  the  logarithmic  versions  of  PLCC  and  PLVC  also 
produce  very  similar  luminance  error. 

The  ANOVA  on  u  V  error  showed  that  the  main  effects 
of  characterization  and  resolution,  plus  their  interaction,  arc 
significant  and  account  for  9.0.  89.9.  and  0.3%  of  the 
variance,  respectively.  Neutral-point  characterization  y  ields 
lower  error  than  conventional  characterization  at  all  resolu¬ 
tions  except  2  and  3.  For  conventional  characterization, 
error  reaches  a  minimum  at  resolution  4  and  then  increases. 


stabilizing  somewhere  around  resolution  9.  although  the 
changes  beyond  resolution  5  are  trivial.  For  neutral-point 
characterization,  which  shows  the  same  pattern,  the  corre¬ 
sponding  numbers  are  7.  9,  and  4. 

Our  conclusions  about  PLGVC  are  the  same  as  for 
PLGCC.  Both  types  of  error  seem  to  stabilize  around  reso¬ 
lution  9.  so  there  is  no  advantage  to  using  higher  resolu¬ 
tions.  The  CIELUV  and  CIELAB  values  and  their  95th 
percentiles  in  Table  I  suggest  that  resolutions  4  or  5  are  the 
best  choices;  our  tests  with  a  different  monitor  and  colori¬ 
metric  instrument  indicate  that  this  idea  is  gcneralizable.  We 
recommend  9,  however,  for  those  who  prefer  a  conservative 
alternative.  Based  on  the  CIELUV  and  CIELAB  values,  wc 
also  conclude  that  neutral-point  characterization  is  equal  or 
superior  to  conventional  characterization  at  resolutions  of  4 
or  more. 

Summary  and  Discussion 

For  GOG,  resolutions  beyond  3  did  not  produce  substan¬ 
tial  reductions  in  error.  This  finding  agrees  w  ith  Bcms  el  a/.3 
(see  their  Table  II).  who  found  no  differences  among  reso¬ 
lutions  5.  16.  and  180.  For  PLCC  and  PLVC.  there  were  no 
worthwhile  error  reductions  beyond  resolution  13.  For 
PLGCC  and  PLGVC.  resolutions  4  and  5  produced  the  best 
results,  although  the  inconsistencies  at  the  lower  resolutions 
might  justifiably  cause  some  users  to  prefer  resolution  9. 

For  all  resolutions  that  produced  acceptably  small  errors, 
neutral-point  characterization  produced  less  error  than  con¬ 
ventional  characterization.  This  finding  contradicts  the  re¬ 
sults  for  GOG  reported  by  Bcms  el  at.*  (sec  their  Table  II). 
which  showed  no  effect  of  characterization  method.  An¬ 
other  noteworthy  point  is  that  the  magnitudes  of  the  errors 
we  obtained  for  GOG  and  PLCC  are  larger  than  those 
reported  by  Bems  et  al.  (sec  their  Table  III).  Wc  discuss 
these  differences  after  reporting  the  results  for  Experiment  3. 

EXPERIMENT  2:  BLACK  CORRECTION 

Our  second  experiment  examined  the  effects  of  using  black 
correction.  The  black  light,  as  measured  by  the  LMT.  was 
0.06  cd/nr.  u'  =  0.0280.  and  v'  =  0.4529. 

The  procedure  for  black  correction  is  fairly  straightfor¬ 
ward  and  identical  to  correcting  for  ambient  illumination: 
(I)  Zero  the  measuring  target's  DAC  values  while  leaving 
the  background  color  unchanged;  (2)  Measure  the  target's 
tristimulus  values:  (3)  Subtract  these  trislimulus  values  from 
each  measurement  in  the  characterization  dataset;  and  (4) 
When  predicting  the  tristimulus  values  produced  by  a  given 
DAC-valuc  triplet,  do  the  calculation  in  the  usual  way  and 
then  add  the  tristimulus  values  obtained  for  the  target  at  zero 
DAC  value.  Step  3  must  be  performed  for  both  the  conven¬ 
tional  and  neutral-point  characterization  datasets,  when  us¬ 
ing  PLVC  or  PLGVC  with  neutral-point  characterization. 
For  predictive  methods  that  assume  constant  channel  chro- 
maticity.  it  is  important  to  ensure  that  the  assumed  channel 
chromalicities  are  corrected,  as  well  as  the  characterization 
dataset. 
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TABLE  II.  Mean  errors  from  Experiment  2. 


Conventional  characterization  Neutral-point  characterization 


Method 

Rea 

Black 

|%Y| 

uV 

AC 

95ptl 

AC 

95ptl 

l%Y| 

uV 

AC 

95ptl 

AC 

95ptl 

GOG 

4 

N 

1.20 

0.00147 

1.19 

2.68 

0.96 

2.08 

0.66 

0.00086 

0.80 

1.44 

0.63 

1.23 

Y 

0.90 

0.00084 

0.76 

1.33 

0.60 

0.94 

0.52 

0.00070 

0.69 

1.20 

0.52 

0.96 

7 

N 

1.21 

0.00144 

1.17 

2.51 

0.94 

1.96 

0.64 

0.00084 

0.78 

1.39 

0.62 

1.19 

Y 

0.90 

0.00084 

0.77 

1.36 

0.61 

0.96 

0.51 

0.00068 

0.67 

1.20 

0.51 

0.95 

13 

N 

1.20 

0.00142 

1.16 

2.48 

0.93 

1.92 

0.63 

0.00084 

0.78 

1.41 

0.62 

1.18 

Y 

0.90 

0.00082 

0.76 

1.32 

0.60 

0.95 

0.50 

0.00068 

0.67 

1.21 

0.51 

0.94 

PLCC 

7 

N 

1.84 

0.00179 

1.46 

2.85 

1.19 

2.19 

0.87 

0.00131 

1.11 

2.01 

0.87 

1.66 

Y 

1.54 

0.00117 

1.02 

1.91 

0.82 

1.34 

0.78 

0.00116 

1.01 

1.83 

0.76 

1.39 

13 

N 

1.26 

0.00140 

1.15 

2.40 

0.93 

1.87 

0.62 

0.00085 

0.77 

1.43 

0.61 

1.23 

Y 

0.96 

0.00083 

0.77 

1  39 

0.61 

0.99 

0.50 

0.00070 

0.67 

1.22 

0.51 

0.97 

PLGCC 

4 

N 

0.78 

0.00121 

1.00 

2.46 

0.79 

1.90 

1.05 

0.00052 

0.63 

1.14 

0.51 

0.92 

Y 

0.66 

0.00065 

0.65 

1  26 

0.50 

1.01 

0.97 

000046 

0.55 

0.99 

0.45 

0.65 

7 

N 

1.05 

0.00126 

1.06 

2.24 

0.84 

1.70 

0.70 

0.00065 

0.66 

1.23 

0.52 

0.99 

Y 

0.77 

0.00071 

0.71 

1.30 

0.56 

0.98 

0.59 

0.00051 

0.57 

1.05 

0.43 

0.85 

13 

N 

1.11 

0.00131 

1.09 

2.29 

0.87 

1.75 

0.64 

0.00073 

0.70 

1.32 

0.56 

1.05 

Y 

0.82 

0.00076 

0.73 

1.31 

0.57 

0.97 

0.53 

0.00058 

0.61 

1.13 

0.46 

0.90 

PLVC 

7 

N 

1.84 

0.00125 

1.06 

1.97 

0.85 

1.46 

0.80 

0.00104 

0.92 

1.76 

0.73 

1.43 

Y 

1.54 

0.00095 

0.85 

1.47 

0.70 

1.09 

080 

0.00105 

0.95 

1.84 

0.72 

1.43 

13 

N 

1.26 

0.00089 

0.80 

1.52 

0.64 

1.11 

0  56 

0  00059 

0.64 

1.21 

0.52 

1.00 

Y 

0.96 

0.00056 

0.56 

1.10 

0.47 

0.86 

0.59 

0.00060 

0.64 

1.23 

0.48 

096 

PLGVC 

4 

N 

0.78 

0.00065 

0.68 

1.41 

0.55 

1.14 

1.00 

0.00049 

0.66 

1  26 

0.57 

1.20 

Y 

0.66 

0.00057 

0.58 

1.16 

0.48 

1.01 

1.04 

0.00042 

0.57 

1.12 

0.45 

0.86 

7 

N 

1.05 

0.00075 

0.73 

1.42 

0.59 

1.08 

0.65 

0.00046 

0.60 

1.15 

0.49 

0.98 

Y 

0.77 

0.00046 

0.53 

1.06 

0.44 

0.84 

0.66 

0.00046 

0.55 

1.07 

0.41 

0.78 

13 

N 

1.11 

0.00080 

0.75 

1.44 

0.60 

1.07 

0.61 

0.00050 

0.61 

1.16 

0.50 

0.97 

Y 

0.82 

0.00048 

0.52 

1.07 

0.43 

0.82 

0.61 

0.00050 

0.59 

1.12 

0.44 

0.84 

To  ease  the  computational  burden,  we  restricted  attention 
to  resolutions  4,  7.  and  13  for  GOG.  PLGCC.  and  PLGVC. 
and  resolutions  7  and  13  for  PLCC  and  PLVC.’  We  per¬ 
formed  ANOVAs  using  characterization  method,  resolu¬ 
tion.  and  black  correction  (present  vs.  absent)  as  the  inde¬ 
pendent  variables,  and  again  used  an  alpha  of  0.05  for  all 
tests.  The  mean  luminance  and  u  V  errors  arc  summarized 
in  Table  II  (wherein  the  values  for  non-black-corrcctcd 
eases  duplicate  the  corresponding  values  in  Table  I)  and 
illustrated  in  Fig.  2. 

The  ANOVA  results  for  effects  that  were  tested  in  Ex¬ 
periment  I  match  expectations  in  all  important  respects. 
One  new  tiling  we  learned  is  that  the  main  effect  of  black 
light  is  significant  in  all  cases  (emir  was  lower  when  black 
correction  was  used).  The  most  interesting  outcome, 
though,  was  that  the  characterization  x  black-light  interac¬ 
tion  is  significant  in  all  cases,  accounting  for  1.3%  (PLCC) 
to  2().I%(PLGCC)  of  the  variance  in  luminance  error  and 
10.4%  (PLCC)  to  15.2%  (PLGVC)  of  the  variance  in  u’v’ 
error.  The  interaction  is  visible  in  Fig.  2:  Black  correction 
produces  a  much  smaller  emir  reduction  (or  sometimes 
none)  when  neutral-point  characterization  is  used.  Simi¬ 
larly,  neutral-point  characterization  is  less  effective  when 
black  correction  is  used.  These  findings  imply  that  there  is 


'  Given  the  Experiment  I  results  for  PLGCC  anil  PLGVC.  revolution  9 
might  veem  ■  better  choice  than  resolution  7:  however.  Experiment  I 
showed  dial  resolution  1 3  should  give  the  same  results  as  9  for  PLGCC  and 
PLGVC,  so  wc  thought  7  might  be  mote  informative. 


overlap  in  the  error  accounted  for  by  neutral-point  charac¬ 
terization  and  black  correction — an  implication  we  explored 
in  the  next  experiment. 

EXPERIMENT  3:  NUMBER  OF  CHANNELS 

Our  124  test  colors  consisted  of  12  single-channel  colors 
(i.e.,  colors  produced  with  only  one  nonzero  DAC  value), 
48  two-channel  colors,  and  64  three-channel  colors;  thus,  it 
was  weighted  heavily  with  two-  and  especially  three-chan¬ 
nel  colors.  Furthermore,  the  results  from  Experiment  2 
suggest  that  the  (seemingly  small  amount  of)  black  light  had 
a  substantial  effect  on  prediction  accuracy,  and  the  results 
reported  in  Appendix  B  show  that  our  monitor  was  some¬ 
what  sub-additive.  Three  mechanisms  could  cause  these 
facts  to  produce  the  characterization  x  black-light  interac¬ 
tion  we  observed  for  luminance  error  in  Experiment  2,  as 
well  as  the  apparent  superiority  of  neutral-point  character¬ 
ization  over  conventional  characterization  in  Experiment  I: 

1 .  Euch  measurement  for  a  conventional  characterization 
dataset  includes  any  black  light  that  may  be  present. 
Therefore,  if  conventional  characterization  is  used 
without  bluck  correction,  luminance  predictions  for  sin¬ 
gle-channel  colors  should  be  accurate  nonetheless, 
whereas  predictions  for  two-  and  three-channel  colors 
should  ovcr-rcprescnt  the  black  light  by  two  or  three 
times  and.  therefore,  be  less  accurate.  The  situation  for 
neutral-point  characterization  is  the  opposite;  If  black 
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FIG.  2.  Mean  %Y|  and  uV  error  from  Experiment  2  for  each  predictive  method  as  a  function  of  resolution,  characterization 
method,  and  black  correction. 


correction  is  not  performed,  decomposing  the  A'  mea¬ 
surements  distributes  the  black  light  across  the  RiiH 
channels  in  the  characterization  dataset,  so  luminance 
predictions  involving  three  channels  should  ho  more 
accurate  than  predictions  involving  one  or  two.  which 
should  underestimate  the  black  light. 

2.  A  similar  principle  should  apply  if  the  monitor's  addi¬ 
tivity  is  imperfect,  even  if  black  correction  is  used.  In 
this  case,  luminance  predictions  from  conventional 
characterization  should  be  most  accurate  for  single- 
channel  colors,  for  which  additivity  is  irrelevant,  and 
should  become  less  accurate  as  the  number  of  channels 
land  their  DAC  value)  increases  and.  hence,  additivity 
problems  manifest  themselves.  Neutral-point  character¬ 
ization.  on  the  other  hand,  is  derived  from  three-channel 
measurements  and  should,  therefore,  be  most  accurate 
for  three-channel  colors. 

3.  The  luminance  added  by  black  light  is  a  constant: 


hence,  regardless  of  characterization  method,  black- 
light-rclaled  luminance  error  assessed  as  a  percentage 
should  diminish  as  luminance  increases.  In  our  dataset, 
luminance  lends  to  increase  with  the  number  of  chan¬ 
nels  in  use.  and  thus — for  cases  where  black  correction 
was  not  used  percent  luminance  error  should  decrease 
as  the  number  of  channels  increases. 


To  summarize:  For  our  dataset,  mechanisms  I  and  2 
should  cause  luminance  error  to  increase  as  the  number  of 
channels  increases  for  conventional  characterization,  and 
decrease  as  the  number  of  channels  increases  for  neutral- 
point  characterization.  Mechanism  I  should  cause  black 
correction  u*  reduce  error  lor  two-  and  thrcc-channel.  but 
not  single-channel,  predictions  for  conventional  character¬ 
ization  and  be  least  effective  for  three-channel  predictions 
for  neutral-point  characterization.  Mechanism  3  should 
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FIG.  3.  Mean  %Y|  and  u'v'  error  from  Experiment  3  for  each  predictive  method  as  a  function  of  characterization  method, 
black  correction,  and  number  of  channels. 


cause  black  correction  to  also  reduce  differences  caused  by 
the  number  of  active  channels. 

As  for  u'v’  error,  many  complex  mechanisms  occurred  to 
us  that  might  contribute  to  the  interaction.  Mainly,  though, 
we  expected  the  same  mechanisms  that  affect  luminance 
error  to  also  affect  u’v'  error,  because  accurate  chroniaticity 
predictions  often  require  accurate  predictions  of  the  constit¬ 
uent  channel  luminances.  Overall,  then,  it  seemed  that  the 
test  dataset  w-as  weighted  in  favor  of  neutral-point  charac¬ 
terization.  especially  in  the  cases  where  black  correction 
wus  not  performed. 

The  results  in  Table  II  seem  to  lit  our  reasoning  fairly 
well,  so  we  tested  it  by  breaking  down  the  resolution  13  data 
in  Table  II.  according  to  the  number  of  channels  involved. 
Wc  performed  ANOVAs  using  characterization  method, 
black  correction,  and  number  of  channels  as  the  independent 
variables,  and  again  used  an  alpha  of  0.05  for  all  tests.  The 
results  are  summarized  in  Tabic  III  and  illustrated  in  Fig.  3. 


The  ANOVA  results  are  very  similar  for  all  predictive 
methods.  For  luminance  error,  the  main  effect  of  black  light 
and  the  black-light  x  channels,  characterization  x  channels, 
and  characterization  .r  black-light  a  channels  interactions  arc 
significant  in  all  cases;  the  characterization  x  black-light 
interaction  is  also  significant  for  PLVC  and  PLGVC.  The 
characterization  x  channels  interaction  is  consistently  the 
most  important  effect,  accounting  for  72.1%  (GOG)  to 
83.3%  (PLGVC)  of  the  variance.  For  u'v'  error,  all  main 
effects  and  interactions  are  significant  in  all  cases.  For 
GOG.  PLCC.  and  PLGCC.  the  main  effects  of  character¬ 
ization.  black  light,  and  channels  are  the  most  important, 
accounting  for  an  average  of  14.5,  39.5.  and  21.5%  of  the 
variance,  respectively.  For  PLVC  and  PLGVC,  however, 
the  main  effect  of  channels  is  by  far  the  most  important, 
accounting  for  80. 1  and  73.4%.  respectively. 

Examination  of  Table  111  and  Eig.  3  shows  that  luminance 
error  behaves  almost  exactly  as  wc  predicted.  The  one 
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TABLE  III.  Mean  errors  from  Experiment  3. 


1 

Black 

Chn 

Conventional  characterization 

Neutral-point  characterization 

l%Y1 

u'v' 

AS. 

95ptl 

ASu 

95ptl 

|%Y| 

UV 

AS, 

95pll 

AS* 

95ptl 

GOG 

N 

1 

046 

0  00151 

0.66 

1.76 

0.91 

266 

1.58 

0.00151 

0.80 

1.46 

080 

2.26 

2 

0.97 

0.00139 

1.09 

2.67 

080 

1.99 

0.64 

0.00075 

0.76 

1.38 

059 

1.13 

3 

1.52 

000144 

1.31 

2.41 

1.04 

1.75 

0.45 

0.00078 

0.78 

1.41 

060 

1.13 

Y 

1 

0.42 

0.00091 

0.50 

1.26 

0.53 

1.29 

0.92 

000091 

066 

1.11 

0.66 

1.18 

2 

0.66 

0.00074 

0.67 

1.35 

0.50 

0.91 

0.64 

0.00068 

0.67 

1.23 

0.50 

0.93 

3 

1.16 

0.00086 

0.87 

1.32 

069 

0.95 

0.31 

0.00064 

0.66 

1.20 

0.48 

0.87 

PLCC 

N 

1 

0.46 

0.00151 

0.69 

1.97 

0.93 

2.86 

1.52 

0.00151 

080 

1  63 

081 

2.33 

2 

1.03 

0  00136 

1.08 

2.56 

0.80 

1.97 

063 

000078 

0.76 

1.38 

060 

1.11 

3 

1.57 

0.00141 

1.29 

234 

1.02 

1.74 

0.43 

0.00078 

0.77 

1.44 

0.59 

1.21 

Y 

1 

0.50 

0.00091 

0.51 

1.30 

0.54 

1.34 

0  98 

0.00090 

0.68 

1.16 

0  66 

1.23 

2 

0.72 

0.00076 

0.69 

1.38 

0.52 

0.95 

0.65 

0.00071 

0.68 

1.22 

051 

096 

3 

1.22 

000086 

088 

1.40 

0.69 

0.99 

030 

0.00065 

0.66 

1.23 

0.48 

0.90 

PLGCC 

N 

1 

036 

0.00151 

0.60 

1.69 

0.87 

2.61 

1.71 

0.00151 

0.80 

1.42 

0.79 

223 

2 

0.87 

0.00126 

1.01 

2.47 

0.73 

1.81 

071 

0.00065 

0.71 

1.30 

055 

1.02 

3 

1.44 

0.00131 

1.23 

2.22 

0.96 

1.62 

0.39 

000064 
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1.26 
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1.01 

Y 

1 

0.37 

000091 

0.49 
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1.23 

1.01 

0.00091 

0.69 

1.13 
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1  16 

2 

058 

0.00068 

064 

1  30 

0.47 

0.89 

0.71 

0.00059 

0.63 

1.15 

0  48 

0.87 

3 

1  08 

000078 

085 

1.32 

0.66 

0.97 

0.31 

0.00052 

0.58 

1.10 

0.41 

0.78 

PLVC 

N 

1 

0.48 

0.00014 

0  28 

0.59 

0.23 

043 

1.04 

0.00013 

0.59 

1.26 

0.43 

1.02 

2 

1.03 

0.00090 

064 

1  37 

0.51 

0.94 

068 

0.00071 

0.71 

1.23 

0.63 

1.04 

3 

1  57 

000103 

1  01 

1.57 

082 

1.15 

0.39 

0.00059 

0.59 

1.18 

0.46 

0.96 

Y 

1 

0.50 

0.00015 

0.30 

0.59 

0  25 

0.49 

1.07 

0.00015 

0.63 

1.12 

0.43 

0.86 

2 

0.72 

0.00066 

0  55 

1.12 

0  44 

0.81 
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0.00071 

0.69 

1.32 

0.52 

0.99 

3 

1.22 

0.00057 

0.62 

1.12 
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088 

0.38 

000060 

0.60 

1.19 

0.46 

0.95 
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N 

1 

0.36 

0.00014 

0.23 

0.55 

0.19 

040 

1.29 

000014 

0.70 

1.30 

0.50 

1.01 

2 

0.87 

0.00080 

059 

1.30 

048 
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0.75 

0.00063 

0  72 

1.22 

0.65 

1.07 

3 

1.44 

0.00092 

096 

1.52 

0  77 

1.13 

0.38 

0.00046 

0.52 

1.04 

0.39 

0.83 

Y 

1 

0,37 

0.00014 

0.24 

0.55 

0.20 

042 

1.13 

0  00015 
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1.12 

0.44 

0.88 

2 

0.58 

0.00058 

0.50 

1.08 

0.41 

0.75 

0.79 

0  00063 

065 

1.18 

0.50 

0.92 

3 

1.08 

0.00047 

0.59 

1.09 

0.49 

0.87 

0.37 

0.00047 

053 

1  00 

0.40 
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exception  is  the  fact  that,  for  GOG.  PLCC.  and  PLGCC 
using  neutral-point  characterization,  black  correction  pro¬ 
duces  a  greater  improvement  for  three-channel  colors  than 
for  two-channel  colors.  This  exception  suggests  that  there  is 
a  fourth,  as-yet  unidentified  mechanism  influencing  the  re¬ 
sults. 

The  results  for  uV  error  indicate  that,  contrary  to  our 
expectations,  error  in  luminance  prediction  docs  not  neces¬ 
sarily  produce  concomitant  error  in  chromaticity  prediction. 
For  CiOG.  PLCC.  and  PLGCC.  error  for  single-channel 
colors  is  noticeably  greater  than  for  two-  and  three-channel 
colors  when  neutral-point  characterization  is  used  without 
black  correction.  Otherwise,  the  effect  of  the  number  of 
channels  is  modest.  We  do  not  understand  the  latter  out¬ 
come.  but  offer  the  following  explanation  for  the  former: 
When  black  light  is  present  but  unconnected,  decomposing 
the  neutral-point  measurements  places  erroneous  luminance 
values  in  the  characterization  dataset,  due  to  the  use  of 
channel  chromaticilies  in  the  decomposition  calculations 
that  are  correct  only  at  Dm„  (assuming  they  were  measured 
al  Therefore,  for  single-channel  colors,  the  predicted 
chromaticilies  are  usually  wrong,  because  the  assumed 
channel  chromaticilies  are  wrong,  hut  for  two-  and  three- 
channel  colors  the  luminance  and  chromaticity  emirs  offset 
each  other,  yielding  more  accurate  chromaticity  predictions. 

Another  interesting  point  that  is  evident  in  Fig.  3  is  that 
PLVC  and  PLGVC  yield  essentially  the  same  single-chan¬ 


nel  uV  error,  regardless  of  characterization  method  or 
black  correction,  and  this  error  is  substantially  lower  than 
the  other  predictive  methods  in  all  cases.  These  outcomes 
imply  that  the  changes  in  our  measured  channel  chromatici¬ 
lies  were  caused  by  something  more  than  black  light,  and 
the  variable  channel -chromaticity  feature  of  PLVC  and 
PLGVC  docs  a  good  job  of  accounting  for  it.  as  well  as 
black  light,  at  least  for  single-channel  colors.  Although  error 
is  greater  for  two-  and  three-channel  colors,  it  is  often  less 
than  for  the  other  predictive  methods  when  black  correction 
is  used  and.  in  this  case,  the  complications  of  using  neutral- 
point  characterization  in  conjunction  with  PLVC  and 
PLGVC  do  not  appear  to  be  worth  the  effort. 

We  conclude  that  lire  superiority  of  neutral-point  charac¬ 
terization  over  conventional  characterization  that  we  ob¬ 
served  in  Experiment  I  was  an  artifact,  caused  by  the 
presence  of  uncorrccled  black  light,  non-additivity,  and  a 
test  dataset  that  emphasized  multichannel  colors.  Although 
Bcrns  cl  til. '  used  an  equivalent  dataset,  they  stated  that  their 
measuring  geometry  excluded  black  light,  and  evidently 
their  additivity  was  better  than  ours  (although  we  must  note 
that  their  Table  IV  suggests  that  it  was  similar).  If  additivity 
is  good  and  black  light  is  either  absent  or  corrected  for. 
neutral-point  and  conventional  characterization  should  yield 
comparable  results,  regardless  of  the  number  of  channels 
in  use. 

If  the  black-corrected  GOG  and  PLCC  results  for  neutrul- 


Volume  25,  Number  2,  April  2000 


99 


w<fftodRMCrw 


PLCC  13  C 
OOO  7  C 
COO  1)  c 
OOO  4  C 
PLGCC13  C 
PtOCC  7  C 
OOO  4  K 
PICC  13  N 
COO  13  N 
OOO  7  N 
PtOCC  4  C 
PLVC  13  N 
PLOVC  4  C 
PLVC  13  C 
PtOCC  13  N 
PtQVC  4  N 
PtOCC  4  N 

pi  oven  h 
PL CVC  7  C 
PtOCC  7  N 
PLOVC13  C 
PtOVC  7  N 


0  0  02  0.4  00  0  8  10 


AE-ab 


C  *  CwwrJttontl  cft4rictii»aon 
N  ■  Neural  ctunctartrjsxm 

FIG.  4.  Mean  and  95th-percentile  A E^,  for  black-corrected 
cases  from  Experiment  2,  sorted  in  descending  order.  PLCC 
and  PLVC  at  resolution  7  are  omitted. 


point  characterization  in  Table  II  arc  compared  with  those 
reported  by  Bents  el  til ,3  in  their  Tables  II  and  III.  the 
cross-laboratory  differences  arc  much  less  than  the  ones  we 
reported  for  Experiment  I .  Indeed,  the  differences  are  prob¬ 
ably  ns  small  ns  could  be  expected  for  cross-laboratory 
comparisons. 


CONCLUSIONS 


for  a  predictive  method.  This  choice  requires  a  total  of 
only  six  measurements  (three  .V  plus  one  each  of  R,  G. 
and  B  to  establish  the  channels'  chromaticity  coordi¬ 
nates)  or  nine  measurements,  depending  on  whether 
neutral-point  characterization  is  used,  and  yields  very 
good  accuracy. 

4.  If  programming  and  usage  complexity  must  be  mini¬ 
mized.  PLCC  at  resolution  13  is  the  best  choice.  It 
requires  cither  16  or  39  measurements  (depending  on 
whether  neutral-point  characterization  is  used),  is  the 
easiest  predictive  method  to  implement  in  software, 
does  not  require  human  intervention  to  ensure  proper 
convergence,  and  yields  accuracy  comparable  with 
GOG. 

5.  The  best  compromise  between  complexity  and  measur¬ 
ing  time  is  PLGCC  at  resolution  4.  It  is  only  slightly 
more  difficult  to  program  than  PLCC.  requires  either  7 
or  15  measurements  (depending  on  whether  neutral- 
point  characterization  is  used),  and  is  a  bit  more  accu¬ 
rate  than  GOG  or  PLCC. 

6.  If  programming  complexity  is  not  an  issue  and  accu¬ 
racy  must  be  maximized.  PLGVC  using  conventional 
characterization  at  resolution  9  is  the  best  choice.  It 
requires  27  measurements,  should  yield  the  same  accu¬ 
racy  as  resolutions  13  or  higher,  and  should  also  yield 
accuracy  that  is  belter  than  or  comparable  with  PLGVC 
with  neutral-point  characterization,  regardless  of  reso¬ 
lution,  while  avoiding  the  complications  and  added 
measurements  required  to  use  neutral-point  character¬ 
ization  in  conjunction  with  PLGVC.  However,  the  ac¬ 
curacy  advantage  over  the  other  methods  is  realized 
only  for  colors  at  or  near  the  edges  of  the  display’s 
color  gamut. 


To  facilitate  comparisons  of  predictive  accuracy.  Fig.  4 
graphs  mean  A E*ah  in  descending  order  for  the  black-cor¬ 
rected  cases  from  Experiment  2  (minus  PLCC  and  PLVC  at 
resolution  7.  because  these  cases  are  mil  recommended  for 
use).  Tlte  95th  percentiles  are  shown,  too.  Based  on  this 
figure  and  consideration  of  our  other  results,  we  offer  the 
following  recommendations: 

1.  Black  correction  should  be  performed  whenever  there 
is  reason  to  believe  the  area  of  interest  on  the  display  is 
luminous  at  zero  DAC  value.  If  a  colorimetric  instru¬ 
ment  having  sufficient  sensitivity  for  this  purpose  is 
unavailable,  the  use  of  neutral-point  characterization 
provides  much,  but  not  all.  of  the  benefits  of  black 
correction. 

2.  Neutral-point  characterization  is  to  be  preferred  over 
conventional  characterization,  except  for  PLVC  and 
PLGVC  and  in  special  cases  where  most  or  till  of  the 
colors  are  produced  by  single  channels.  Otherwise, 
neutral-point  characterization  yields  equal  or  better  pre¬ 
dictive  accuracy  and  reduces  the  number  of  character¬ 
ization  measurements  needed. 

3.  If  tlte  number  of  characterization  measurements  must 
be  minimized.  GOG  at  resolution  3  is  the  best  choice 


Our  results  indicate  that,  for  most  purposes,  there  is  little 
reason  to  prefer  one  method  over  another  on  the  basis  of 
predictive  accuracy,  given  that  each  is  used  in  accordance 
with  our  recommendations.  Of  particular  interest  is  the  fact 
that  PLVC's  and  Pl.GVC’s  allowance  for  variable  channel 
chromaticity  yields  substantial  advantages  only  lor  single- 
channel  colors.  This  finding  contradicts  our  earlier  compar¬ 
isons'  against  PLCC  and  may  reflect  differences  in  the 
monitors:  therefore,  we  are  uncertain  how  gcncralizablc  this 
finding  is.  We  also  note  that  PLVC  and  PLGVC  may  be 
especially  valuable  for  predicting  the  colorimetry  of  liquid- 
crystal  displays,  which  as  Silverstcin  and  Fiske"  have 
shown,  can  exhibit  substantial  covariance  between  lumi¬ 
nance  and  channel  chromaticity. 

It  is  worthwhile  to  consider  the  adequacy  of  the  predic¬ 
tive  accuracies  we  have  reported.  Stokes.  Fairchild,  and 
Bents'1 2  found  that,  for  sequential  comparisons  of  images  on 
the  same  display,  the  A£*,,  that  yielded  a  50%  discrimina¬ 
tion  threshold  ranged  from  1.43-2.65  and  averaged  2.15. 
All  the  methods  we  studied  can  provide  95th-pcrccntilc 
accuracies  that  are  within  these  limits.  For  simultaneous 
comparisons  of  images  on  adjacent  displays,  on  the  other 
hand — where  the  predictive  errors  can  sum  and  thresholds 
arc  presumably  smaller — adequacy  is  less  clear.  We  can 
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estimate  limiting  thresholds  for  this  case  by  referring  to 
MacAdam's'-'  data:  Brainard14  has  determined  recently  that 
Nutting's  509r  discrimination  thresholds,  expressed  as 
A E^,,  ranged  from  0.30-3.37  and  averaged  1.23.  Therefore, 
none  of  the  methods  we  studied  may  be  fully  adequate  for 
matching  images  across  adjacent  displays.  In  such  cases,  the 
ineasure-und- adjust  algorithm  we  introduced  previously'  may 
be  useful. 

Our  results  suggest  that  the  last  signilicant  emir  source 
remaining  to  be  eliminated  via  modeling  is  non-additivity. 
Although  progress  toward  solving  this  problem  has  been 
reported,*-'^17  a  major  difficulty  is  devising  a  way  to  in¬ 
corporate  channel  interactions  in  the  model  while  preserv¬ 
ing  reversibility  (so  predictions  can  be  made  of  the  DAC 
values  needed  to  produce  a  given  color,  as  well  us  the  color 
produced  by  given  DAC  values)  and  without  making  the 
calculations  forbiddingly  complex.  Our  dataset  is  available 
by  request  to  anyone  who  wishes  to  pursue  this  problem  or 
test  oilier  ways  of  improving  on  the  accuracies  we  have 
reported. 

APPENDIX  A:  Pl.VC  AND  PLGVC 

Our  earlier  article1  explained  PLVC  by  describing  how  it  is 
used  to  calculate  the  DAC  values  that  are  needed  to  produce 
a  desired  luminance  and  pair  of  chromaticity  coordinates. 
Communications  with  other  workers  have  since  indicated 
this  explanation  was  not  sufficiently  clear.  We  attempt  to 
rectify  this  problem  here  by  describing  how  PLVC  is  used 
to  calculate  the  luminance  and  chromaticity  coordinates  that 
arc  produced  by  a  given  set  of  DAC  values.  This  case  is 
much  easier  to  explain  and  understand.  Copies  of  our  PLVC 
and  PLGVC  subroutines,  coded  in  FORTRAN,  arc  also 
available  by  request. 

Assume  that  XLZ-tristimulus  values  have  been  measured 
for  R.  G,  and  B  at  several  DAC  values  each,  with  only  one 
nonzero  DAC  value  at  a  time.  This  yields  a  table  of  mea¬ 
sured  tristimulus  values  for  each  channel.  To  predict  the 
tristimulus  values  that  will  be  produced  for  DH  —  i,  Dc 
0.  D„  -  0,  where  i  >  0.  enter  the  table  for  R  and  interpolate 
between  the  two  measurements  that  lie  immediately  above 
and  immediately  below  i  to  get  the  predicted  XK.  Repeat  this 
procedure  to  get  YR  and  V.t(.  These  predicted  tristimulus 
values  for  R  are  the  final  solution,  because  Da  and  Dtl  are 
zero,  which  we  assume  to  mean  that  the  tristimulus  values 
for  the  G  and  H  channels  are  also  zero.  This  solution 
technique  differs  from  PLCC  in  that  three  independent 
piecewise  linear  interpolations  are  used  to  obtain  the  three 
tristimulus  values.  Thus,  variations  in  the  measured  R  chro¬ 
maticity  coordinates  as  the  R  DAC  value  changes  are  ac¬ 
counted  for  in  the  resulting  predictions.  With  PLCC.  on  the 
other  hand,  interpolation  is  used  only  to  obtain  the  predicted 
y-lristimulus  value.  The  predicted  X  and  Z  are  obtained 
from  the  predicted  Y  plus  a  single  pair  of  measured  chro¬ 
maticity  coordinates  for  the  R  channel,  which  are  assumed 
to  be  constant. 

A  slightly  more  complicated  example  involves  DK  —  i, 
Da  =  j,  D„  k,  where  [/,  j,  k)  >  0.  In  this  case,  the 


procedure  described  above  must  be  performed  for  each 
channel,  resulting  in  a  total  of  nine  piecewise  linear  inter¬ 
polations.  Each  interpolation  yields  a  predicted  tristimulus 
value.  The  final  solution  is  obtained  by  summing  these 
tristimulus  values  across  each  channel;  that  is.  computing  X 
—  XR  +  X(;  +  X„.  etc.  PLGVC  works  the  same  way  as 
PLVC.  but  interpolation  is  performed  for  logtX*)  as  a 
function  of  logfD*),  etc.  Thus,  PLVC  and  PLGVC  perform 
trilincar  interpolation  based  solely  on  measurements  of  sin¬ 
gle  channels,  and  arc  similar  with  (but  simpler  than)  trilincar 
interpolation  metluxls  that  include  multichannel  measure¬ 
ments  and  have  become  popular  for  predicting  the  output  of 
color  printers. 

In  our  implementation  of  PLVC  (und  PLCC.  for  that 
matter),  predictions  for  DAC  values  that  lie  below'  the 
smallest  measured  values  in  the  tables  are  obtained  by 
assuming  that  the  tristimulus  values  become  zero  at  zero 
DAC  value.  That  is.  we  interpolate  between  zero  and  the 
smallest  measured  values.  We  follow  this  procedure  be¬ 
cause  it  is  difficult  ordinarily  to  obtain  reliable  colorimetric 
measurements  at  DR  —  Dc  -  D„  =  0.  especially  because 
we  do  our  best  to  assure  that  the  monitor  does  not  emit  at 
this  setting. 

A  PLVC  algorithm  to  predict  tristimulus  values,  given 
DAC  values,  is  only  slightly  more  difficult  to  implement  in 
software  than  PLCC.  The  opposite  operation — that  is,  pre¬ 
dicting  the  DAC  values  that  yield  a  given  set  of  tristimulus 
values — is  more  troublesome.  One  method  is  iterative 
search;  Compute  XYZ  for  every  possible  combination  of  DH, 
/•>,,.  and  l)„  and  select  the  DK,  Da.  Dt,  triplet  that  comes 
closest  to  yielding  the  desired  XYZ  (using  w  hatever  measure 
of  color  error  one  prefers).  This  method  is  inefficient  and. 
depending  on  the  application  and  the  number  of  bits  in  one’s 
graphics  DACs,  can  be  unucceptably  slow.  Appendix  A  in 
our  previous  article  describes  a  far  more  efficient  algorithm 
(hat  we  have  used  routinely  for  many  years. 

Our  PLVC  algorithm  does  not  yield  valid  solutions  for 
PLGVC.  however.  That  is.  one  cannot  simply  log-transform 
the  DAC  values  and  measured  XYZ,  use  our  PLVC  algo¬ 
rithm  to  find  log(/J*.  D(;,  D„).  und  then  transform  back.  The 
problem  is  that  a  closed-form  solution  to  the  equation  our 
PLVC  algorithm  needs  in  order  to  perform  its  iterations 
does  not  exist  in  the  log-log  domain. 

The  best  solution  to  this  problem  we  have  devised  so  far 
takes  advantage  of  the  facts  that  the  correct  solution  sub- 
space  for  PLVC  must  also  be  the  correct  subspace  for 
PLGVC.  and  furthermore  the  values  of  I)R,  l)G,  and  D„  that 
constitute  the  PLGVC  solution  must  be  greater  than  or  equal 
to  the  values  for  the  corresponding  PLVC  solution.  There¬ 
fore.  we  compute  a  solution  using  PLVC.  thereby  identify¬ 
ing  the  correct  subspace  for  PLGVC,  and  then  perform  an 
iterative  search  of  all  DH.  Da,  D„  triplets  lying  between  the 
PLVC  solution  and  the  top  of  the  solution  subspace.  Thus, 
we  obtain  solutions  using  iterative  search,  but  reduce  the 
required  number  of  iterations  greatly  by  restricting  the 
search  area.  With  this  approach,  the  required  number  of 
iterations  decreases  as  the  resolution  of  the  characterization 
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dataset  increases,  because  the  si/e  of  the  solution  subspaces 

decreases. 

APPENDIX  B:  CORRECTED  NEI  TR  VI -POINT 
CHARACTERIZATION 

There  is  an  oversight  in  the  description  given  by  Bern*  ei 
al/  (and  the  CIE4)  for  implementing  neutral-point  charac¬ 
terization.  One  of  the  equations  used  by  Bents  el  al.  to 
predict  the  CIE  1931  XFZ-lristimulus  values  that  is  pro¬ 
duced  by  specified  DAC  values  is 


X 

Xc.jK",  Xkji*.. 

Y 

= 

Y MJXmrn  Y Y 

n- 

Z 

Zgjk.  Y.njK_ 

L*H 

where  .  etc.  arc  the  XKZ-tristimulus  values  obtained 
from  independent  measurements  of  the  R.  G.  and  B  channels 
at  the  maximum-possible  DAC  value  ( D „*„)  and  f£~\ 
and  arc  RGB  monitor  tristimulus  values,  which 
represent  the  outputs  of  the  RGB  channels  and  arc  normal¬ 
ized  with  respect  to  TKU^.  Tau^.  and  respectively. 

For  purposes  of  the  present  discussion,  it  may  he  convenient 
to  think  of  /#-• ,  etc.  as  normalized  luminances  produced  by 
the  R,  G,  and  B  channels,  but  they  can  be  calculated  using 
any  units  that  are  proportional  with  the  radiances  produced 
by  die  RGB  channels. 

When  GOG  is  used  to  predict  the  XTZ-tristimulus  values 
that  me  produced  by  specified  DAC  values,  values  of 

and  are  obtained  by  evaluating  the  GOG  power 
functions  for  those  DAC  values  and  then  Eq.  ( 1 )  is  used  to 
compute  XYZ.  When  GOG  is  used  to  predict  the  DAC 
values  that  yield  desired  XT/.,  the  inverse  of  Eq.  ( 1 )  is  used 
to  compute  the  values  of  1%““.  and  that  yield 
XYZ.  and  inverse  GOG  power  functions  arc  then  used  to 
compute  the  DAC  values  that  yield  etc. 

If  channel  additivity  is  perfect,  evaluating  the  inverse  of 
Equation  I  for  yields  t%“m  =  ~  =  I. 

Typically,  though,  channel  additivity  is  imperfect,  and  these 
values  differ  from  unity.  In  such  cases,  if  neutral-point 
characterization  is  used.  Bents  el  al.  and  the  CIE  recom¬ 
mend  renormalizing  all  values  of  Vjf-,  /§■*",  and  that 
are  obtained  via  the  inverse  of  Eq.  (1)  with  respect  to  the 
values  obtained  for  Renormalization  yields 


So.  when  neutral-point  characterization  is  used  with 
GOG.  the  measurements  of  N  are  decomposed  using  the 
inverse  of  Eq.  (1).  the  results  are  renormalized  if  necessary 


using  Eq.  (2).  and  coefficients  arc  fit  to  the  GOG  power 
functions  to  predict  the  resulting  monitor  tristimulus  values. 
The  outputs  of  the  power  functions  are  used  subsequently  in 
conjunction  with  the  3  X  3  matrix  in  Eq.  (1)  to  predict  XYZ. 
for  specified  DAC  values. 

This  procedure  produces  an  inconsistency,  however, 
when  renormalization  is  used.  The  inconsistency  can  be 
illustrated  by  substituting  the  renormalized  monitor  tris- 
tintulus  values  into  Eq.  (I).  expanding  the  3  X  3  matrix,  and 
rearranging  terms  to  yield 


X 

A*  AC  XH 

>’*  yc  y* 
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where  xH,  etc.,  are  the  chromalicity  coordinates  obtained 
from  the  individual  measurements  of  R.  G.  and  R.  Thus,  it 
can  be  seen  that  the  use  of  renormalized  monitor  tristimulus 
values  with  Eq.  (1)  implicitly  multiplies  monitor  tristimulus 
values  that  have  been  normalized  with  respect  to  Wnat  by 
Yft,,  ,  Y(;d  .  and  Yan  .  That  is.  monitor  tristimulus 
values  dial  have  been  normalized  w  ith  respect  to  the  peak 
white  arc  multiplied  by  the  peak  luminances  obtained  from 
measuring  the  channels  individually,  rather  than  the  lumi¬ 
nances  the  tristimulus  values  were  actually  normalized 
against  (i.c.,  the  peak  channel  luminances  implied  by  the 
peak  white). 

This  inconsistency  adds  unnecessary  error  to  die  GOG 
predictions.  For  example,  if  Eq.  (3)  is  used  to  predict  the 
XFZ-tristimulus  values  produced  by  ~  ~ 

—  1.  the  resulting  XYZ  arc  the  sums  of  the  values  obtained 
from  the  individual  measurements  of  R,  G,  and  B  at 
rather  than  the  values  obtained  for  If  XYZ  arc  calcu¬ 

lated  for  fjy-*,  and  corresponding  with  other  N 
measurements,  the  results  do  not  match  the  values  that  were 
actually  obtained,  cither. 

One  way  to  correct  the  problem  is  to  replace  Eq.  (I )  with 
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whenever  channel  additivity  is  imperfect,  where 
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and  Xu.  .  etc.  are  the  XFZ-tri  stimulus  values  of 
Alternatively,  replace  Eq.  (1)  with 
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TABLE  IV.  Mean  errors  for  Equations  3  vs.  4  at  resolution  13. 


Method 

Eq. 

IttYj 

uV 

AC 

AC 

GOG 

3 

1.41 

0.00119 

0.83 

1.01 

4 

0.63 

0.00084 

0.62 

0.78 

PLCC 

3 

1.44 

0.00119 

083 

1.01 

4 

0.62 

0.0008S 

0.61 

0.77 

PLGCC 

3 

1.31 

0.00109 

0.77 

0.95 

4 

0.64 

0.00073 

0.56 
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The  choice  depends  on  whether  one  finds  it  easier  to  com¬ 
pute  the  3X3  matrix  in  Eq.  (4)  or  switch  between  normal¬ 
ized  and  unnormalized  monitor  tristimulus  values,  as  re¬ 
quired  by  Eq.  (6).  For  methods  like  PLCC  and  PLGCC. 
which  do  not  require  the  use  of  monitor  tristimulus  values. 
Eq.  (6)  can  be  replaced  with 
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which  is  definitely  easier  to  work  with. 

We  compared  Eqs.  (3)  and  (4)  by  performing  calculations 
both  ways  for  GOG.  PLCC.  and  PLGCC  at  resolution  1 3  for 
all  five  replications.  The  results  arc  summarized  in  Table 
IV.  The  error  reductions  due  to  Eq.  (4)  are  similar  for  all 
three  predictive  methods.  Ibe  values  of  /*"*■.  and 
that  were  obtained  for  lVmiU  averaged  0.99450. 
0.97659.  anil  0.98378.  respectively,  across  the  live  neutral- 
point  characterization  datasets.  Therefore,  it  seems  that  our 
correction  yields  useful  error  reductions,  even  for  seemingly 
small  departures  from  additivity. 

APPENDIX  C:  C(Xi  AND  CHARACTERIZATION 
MEASUREMENTS  AT 

There  is  a  problem  with  using  characterization  measure¬ 
ments  performed  at  in  conjunction  with  GOG.  The 
GOG  power  function  for  predicting  the  monitor  tristimulus 
value  for  a  given  channel  has  the  form 

T-Tmu(k.+  kiD/DmJ\  (8) 

where  D  is  the  selected  DAC  value  for  a  given  color 
channel:  Dma,  is  the  maximum  possible  DAC  value:  Tma,  is 
the  channel's  monitor  tristimulus  value  at  k„.  ks,  and 
y  arc  coefficients  representing  the  channel’s  offset,  gain, 
and  gamma,  respectively:  and  T  is  the  predicted  monitor 
tristimulus  value.  Bems  el  at showed  that  k„  +  kt  =  1; 
therefore,  by  substitution.  Eq.  (1)  can  become 


U I  -k.  +  k^/D^  (9) 

or 

T  =  T^(k„  +  ( I  -  A.)  D/D^)\  (10) 

Either  way,  only  two  free  coefficients  appear,  so  only  two 
measurements  arc  needed  to  uniquely  determine  the  values 
of  these  coefficients  for  each  channel.  To  the  best  of  our 
knowledge,  no  closed-form  solution  is  possible,  so  nonlin¬ 
ear  regression  (i.e.,  search  techniques)  is  the  only  available 
means  of  obtaining  solutions.  We  used  the  Marquardt 
method,  as  implemented  in  the  Statistical  Analysis  System's 
NUN  procedure,  for  all  GOG  tests  reported  in  this  article. 

By  definition,  a  measurement  at  Dmax  yields  Tmiu:  there¬ 
fore.  setting  D  —  Omux  in  Eq.  (9)  or  (10)  yields  the  identity 
=  Tmu.  If  one  attempts  to  solve  for  the  coefficients 
using  only  two  measurements,  one  of  which  w  as  obtained  at 
Dmar,  the  measurement,  therefore,  places  no  con¬ 
straints  on  the  coefficients.  One  is  left  with  two  free  coef¬ 
ficients  to  fit  the  one  remaining  measurement,  and  an  infinite 
number  of  solutions  becomes  available.  If  three  or  more 
measurements  are  used  instead,  only  one  solution  is  possi¬ 
ble.  but  it  is  the  same  one  that  would  have  been  obtained  if 
the  D„U1,  measurement  had  been  omitted.*  In  other  words, 
measurements  at  /)w  contribute  nothing  to  solutions  to 
Eqs.  (9)  and  <  10). 

The  same  problem  can  arise  in  a  more  complicated  form 
with  Eq.  (8).  which  demands  a  minimum  of  three  measure¬ 
ments  to  obtain  u  unique  solution.  As  the  regression  pro¬ 
gram  converges  toward  values  of  k„  and  kg  that  sum  to 
unity,  a  Dnuil  measurement  contributes  progressively  less  to 
the  solution,  and  so  the  estimates  of  k„,  kt.  and  y  become 
unstable.  If  more  than  three  measurements  arc  used,  a  D,m;t 
measurement  does  not  interfere  with  convergence,  but  this 
is  because  it  contributes  to  the  solution  only  in  the  early 
stages. 
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